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The general expression for the magnetoresistance (MR) due to the Lorentz force is derived by 
using the Fermi liquid transport theory based on the Kubo formula. The obtained gauge-invariant 
expression is exact for any strength of the interaction, as for the most singular term with respect 
to l/7fe (7fe being the quasiparticle damping rate). By virtue of the exactness, the conserving laws 
are satisfied rigorously in the present expression, which is indispensable for avoiding unphysical 
solutions. Based on the derived expression, we can calculate the MR within the framework of the 
Baym-Kadanofi' type conserving approximation, by including all the vertex corrections required by 
the Ward identity. The present expression is significant especially for strongly correlated systems 
because the current vertex corrections will be much important. On the other hand, if we drop all 
the vertex corrections in the formula, we get the MR of the relaxation time approximation (RTA), 
which is commonly used because of the simplicity. However, the RTA is dangerous because it 
may give unphysical results owing to the lack of conserving laws. In conclusion, the present work 
enables us to study the MR with satisfying the conserving laws which is highly demanded in strongly 
correlated electrons, such as high-Tc superconductors, or ganic metals, and h eavy Fermion systems. 
In Appendix D, we reply to the comment by O. Narikiyo |cond-mat /0006028 1 . (Note that Appendix 
D exists only in the e-preprint version.) 

PACS numbers: 72.10.Bg, 74.25.Fy, 74.70.-h 



I. INTRODUCTION 



The transport phenomena under the magnetic field are very important and interesting field in condensed matter 
physics. Especially, the magnetoresistance (MR) is of current interest: It has been attracting great attentions of many 
researchers because it gives us rich information on the electronic structure of the system. There are several possible 
mechanisms of the MR, for example, (i) the Lorentz force on conduction electrons (orbital effect), (ii) the double- 
exchange mechanism (spin effect), and (iii) the Anderson localization mechanism in dirty metals, or the Kondo effect 
caused by the magnetic impurities. It is well-known that the huge negative MR due to (ii) is realized in Mn-oxides 
at the ferromagnetic transition temperature, i.e., the colossal magnetoresistance. The negative MR due to (iii) is 
observed in dirty semiconductors. 

On the other hand, in usual pure paramagnetic metals, the positive MR due to (i) is the most dominant, which is 
called the orbital MR. There is a long- history on the study of the orbital MR based on the Boltzmann transport theory 
. In many theoretical studies, the collision term in the Boltzmann equation is simply approximated by introducing 
the phenomenological relaxation time (rk), which is called the relaxation time approximation (RTA). There, the 
correlation effect between particle-hole pair is totally dropped. Although RTA is a quasiclassical approximation, 
various important and useful results have been obtained. For example, the MR as a function of the angle between B 
and / in various metals was reproduced within the RTA [Q. 

Especially, the Kohler's rule is one of the most essential relations derived by the RTA: According to this rule, the 
conductivity ct^j^ in the magnetic field B is represented in the functional form a^^ — r • (tB) , if the anisotropy of 
Tk is not so large. According to the Onsager's reciprocity theorem, Ff^^{x) is an even- function of x as for the diagonal 
conductivity {ji = v), and it is an odd- function as for the off-diagonal conductivity {fj, = x^u = y). 

By definition, the resistance without the magnetic filed (po), the Hall coefficient (-Rh), and the magnetoresistance 
(Ap/po), are given by 

Po = l/o-Q 

i?H = {Aa,y/B)/al (1) 
Ap/po = -Aa^x/a-Q - {Aaxy/crof, 
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where (Jq is the diagonal conductivity without B. Aaxy and Aaxx are components proportional to B and B^, 
respectively. We call Aaxx the magnetoconductivity (MC). As a result, the Kohler's rule predicts the relations 
po oc T, i?H oc r°, and Ap/pQ oc t^-B^. These relations are well satisfied in various ordinary metals, which may be 
interpreted as nearly free electron systems. 

Interestingly, however, the Kohler's rule is strongly violated in several strongly correlated metals according to recent 
experimental studies. For example, i?H in the normal state in high-Tc cuprates is approximately proportional to T~^, 
and Ap/ Po oc ^ ' These seemingly non-Fermi liquid behaviors of the transport phenomena in high-Tc cuprates 
have been studied as a central issue both theoretically and experimentally because they should reflect the fundamental 
electronic property of high-Tc cuprates. 

Up to now, several authors have studied the transport phenomena in high-Tc cuprates based on the RTA, and claimed 
that the strong anisotropy of Tk can explain the observed violation of the Kohler's rule. However, its effect is too 
insufficient to reproduce the prominent non-Fermi liquid behaviors observed for wider range of temperatures. On the 
other hand, we recently studied the i?H in high-Tc cuprates based on the general formula for Aaxy derived by Ref . , 
which is formally exact of order r^. Then, we find out that the vertex corrections for the current, which are dropped 
in the RTA, cause the Curie- Weiss like behavior of i?H in the presence of AF fluctuations 

The aim of this paper is to derive the general expression for the MC from the Kubo formula. The obtained expression 
for Aaxx, which is formally exact in order t^, contains all the vertex corrections which ensure the conserving laws. 
In the system with strong correlations, the role of the vertex corrections will qualitatively change the behavior of the 
MR, as it does for the Hall effect. In this respect, the RTA is unreliable because all the current vertex corrections are 
neglected there. Based on the derived expression, we will study the MR in high-Tc cuprates and discuss the violation 
of Kohler's rule in later publications ||]. 

In the absence of the magnetic fleld, Eliashberg derived in 1961 the expression for the conductivity of the Fermi 
liquid system from the Kubo formula H : 



VB JfS M 



(2) 



which is exact up to 0(t). J^^dSk means the two-dimensional integration on the Fermi surfaces, and vb — i^ir)^ 
in the cubic lattice. Tk = l/(— 2ImSk(+iO)) > 0, where Sk(^) is the self-energy. Wka; = {d/dkx){e^ + ReSk(O)) is 
the quasiparticle velocity without the renormalization factor Zk- Moreover, Ji^x has the four^oint vertex correction 
T22, which is not given by the self-energy through the Ward identity. (Jka; is given by eq.([70|) in this paper.) This 
non-trivial current vertex correction from is indispensable to avoid unphysical results. For instance, Yamada 
and Yosida proved that axx given by eq.(^) diverges unless any Umklapp scattering processes exist, reflecting the 
momentum conservation law [ p^ . This physically reasonable result cannot be reproduced if we drop the vertex 
correction from 722- 

The exact expression for the Hall conductivity of order was given by Kohno- Yamada as follows: |p[ 
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(2rk)^ (3) 



where the uniform magnetic fleld is along the z-axis. Here the carrier of an electron is — e (e > 0). Note that if we 
replace Jk^ with t^k^ in eqs. and (||), we get the results of the RTA. 

As for Aaxy, the vertex corrections in Jk^ also play important roles in strongly correlated systems. In fact, according 
to recent theoretical studies, the vertex corrections in J\^x give rise to the strong enhancement of i?H in nearly AF 
Fermi liquids, like in high-Tc superconductors or in k- (BEDT-TTF) organic superconductors . Noteworthily, 
becomes negative in electron-doped compounds irrespective of the fact that its Fermi surface is hole-like everywhere 
|ri|. This long-standing mystery, which cannot be explained within the RTA, is naturally reproduced in Refs. and 
\% due to the fact that Jk is no more parallel to i7k- 

The purpose of this article is to derive the conserving expression for the magnetoconductivity Aaxx which is exact 
in order . We notice that the magnetoconductivity within the RTA, Aa^]:^, is given by 



Aa 



RTA 

XX 



-B' 
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Avb 



I ^ I (i/k X e,) • V (2Tkt;kx)|' (2rk), 



(4) 



where Cz is the unit vector along the magnetic fleld 0|. Note that the expression (Q) has the contribution from the 
fc-derivative of Tk, while axx and Aaxy given by eqs. (||) and (||) do not have it. Zitr^J^ given by eq. ^ will be 
unreliable in correlated electron systems because conserving laws are violated in the sense of Baym and Kadanoff [ p2[ . 
Whereas, the derived formula for Aaxx in this paper, eq.([77|) or eq.([78|) in §V, is conserving and exact up to Oij^. 
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In this sense, our formula for A(Jxx is indispensable for the study of the strongly correlated Fermi liquid. Note that 
our formula coincides with eq. if all the vertex corrections arising from are neglected. 

The contents of this paper are the following: In §11, we study the tight-binding model with Coulomb interaction, 
where the vector potential is included as the Peierls phase in the hopping integrals. In §111, the expression for the 
magnetoconductivity (MC) is obtained without vertex corrections, for an instructive purpose. It coincides with eq.(^) 
at lower temperatures. In §IV and V, we derive the general expression for the MC by taking all the vertex corrections 
into account, which is the main part of this work. In the former section, we derive all the vertexes corrections by 
taking account of the Ward identities for multi-point vertices seriously. In the later section, we perform the analytic 
continuation and derive the exact expression for the MC of order t^. It is given by eq.j?^) or eq.(|7^). Based on the 
obtained formula, we present several discussions on the vertex corrections, and give some remarks in §VI. 



II. KUBO FORMULA FOR THE MAGNETOCONDUCTIVITY 
A. Tight Binding Model in a Magnetic Field 

The tight-binding model without the magnetic field is given by 

Hb=Q = -^3=0 + -^int 

where t^j is the hopping integral and is the spectrum of the non- interacting particles, which is given by the Fourier 

transformation of y cl^ i'^ka-) creation operator of the electron at (with momentum k) with the spin a. 

Hereafter, we may write the vectors and fc as r and k for simplicity. i?int represents the two-body interaction term. 
In this manuscript, we assume the density-density interactions, like the on-site Coulomb interaction. Note that we 
exclude the processes which would change the form of the current operator, for example, the pair-hopping processes 

Now, we introduce the magnetic field through the vector potential Ai on i-site. We do not consider the Zeeman 
term in this article because we focus only on the orbital effect, i.e., the Lorentz force. If the magnetic field is almost 
uniform, Ai is included in the Hamiltonian through the hopping integral as the Peierls phase factor |p^ : 

= ■ exp( -ie{A, +Aj)- (f, - r,)/2 ), (6) 

Here and hereafter, — e (e > 0) is the charge of an electron, and we put c — h — 1. 
From now on, we introduce the external vector potential as 

and extract the coefficient oi B-B' — {iq x A)-{iq' x A') of the conductivity. In the final stage, we put A ^ A' and 
take the uniform limit q = q' — 0. This procedure was originally developed in studying the Hall conductivity in dirty 
metals by Fukuyama et.al. In this article, we take the z-axis to be the direction of the magnetic field. 

Now we expand the Hamiltonian with respect to A and A' . For the purpose of this paper, we need the terms 
proportional to A, A' and AA' . They are given by 

Hb = Hb=o + eA„j„(-g) + eA'^^ ■ jp{-q') + e^A^A'p ■ jap{~q - q') + ■ ■ ■ , (8) 

which is exact up to 0{q, q') and 0{qq'). Here and hereafter, we promise that the summation with respect to a or /3 
is taken implicitly. In cq.(||), jaip), japip) sue defined as 

ja{p) = '^'"lik)cl_p/^Ck+p/2, (9) 
k 

iafi{p) = E ^'a/3(fc)4-p/2Cfe+p/2, (10) 
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where = dae\ and v^p{k) = da/se^, respectively. (Hereafter, da = d/dka and dap = d"^ /dkadkp, etc.) The 

equation (||) is derived in Appendix A. 

Next, we examine the current operator in the magnetic field. It is given by — j + eA/m for a electron in vacuum 
as is well known. However, it is no more valid for a tight-binding model. By using the kinetic equation for an electron 
at f* in the interacting representation, the current operator of the system, j^{r), is given by 

J^(r)=i[F°,fp(f)], (11) 

where p^r) is the density operator of the system |jl^. Based on this definition, j^{p) in the weak magnetic field is 
derived in Appendix A. According to the result, 

j!^{p^O)=jAO) + eAa-j^a{-q) + eA'p-j^a{-q') + e^AaAyj^ap{-q-q') + (12) 

k 

where v1^p{k) = d^a/s^k- Equation ([l^ ) is exact up to 0{q,q') and 0{qq'). In eqs.(||) and (|l^), we dropped the terms 
proportional to A^, A'^ because they are not required for the purpose of this paper. In the same way, we find 

j^iq + q') - + q) + • j^„(g') + eA'^ ■ j^p{q) + e^AaA'^ ■ j^a^(O) + • ■ • , (13) 
up to 0{A,A') and 0{AA'). 

B. Kubo Formula for the Conductivity 

According to Kubo formula jl^ , the conductivity in the magnetic filed is given by 

(7^^(q + q', OJ)^^ [$^^(g + q',0J + iO) - $^^(q + q' , iO)] , (14) 

<p^Aq + q',^x)^ I dTe-"^^(r,j^(g + g',o)j,^(o,T))B, (15) 







where (3 = l/T and (jJ\ = i7rT-2A (A being the integer), respectively. Here T^. is a r ordering operator. $^^((7+17', cj+iO) 
in eq.(|l5|) is given by the analytic continuation of $pi/(g + q\i^x) in eq.(|lj) with respect to Im^A > 

The diagonal conductivity cr^^ does not contain the term proportional to B^n+i because of the Onsager's reciprocity 



theorem. The diagonal conductivity given by cq.(15) can be expanded with respect to q,q' and A, A' as 



cy^.^.{q + 0) = + 2C^;^^''' • (iq^Aa ■ iq',,A'f,) + (16) 

The coefficient C^^^'^p should satisfy the relation C^P'^p — C^^f^'^^ ■ ^zpa^zp'p because of the requirement from the 

gauge-invariance. We caU Z\ct^^ = C^^''^''' • {iqpAa ■ Wp'^p) = C^p'-^y ■ B"^ the magnetoconductivity (MC). The aim of 

this paper is to derive the expression for C^I^'^p by using the Fermi liquid theory. 

In order to derive C^^'^'^ in eq.(p!6[), we study the term proportional to e^AaA'^ in eq.(|l5[), ^^p}{q + q',uj\): It is 
given by 

/9 



$(,iH'? + 9',c^A)= dTe-^--[{TrJp{q + q',0)j,ap{~q~q',T)) 



+ / dT'{TrJf,{q + q',0)j,a{-q,r)jp{~q',T')) 







+ / dT'((a,g)^(/3,g')) 







rP 

+ I dT'{T^j^{q + q',0)j,{0,T)japi-q-q',T')) 




+ dT'dT"{Trjp{q + q',0)MO,r)ja{-q,T')ji,{-q,T"))} + (17) 
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FIG. 1. (i) The Bethe-Salpeter integral equation for the four point vertex r(fc_; k+\k'^\k'_). is the irreducible vertex, (ii) 
Integral equations for A^, Aq and A/j. (iii) The integral equation for A^/j. 



In later sections, we take the derivative of ^\iv{q + q',(jJ\) with respect to qp and q'p,, and perform the ana- 
lytic continuation with respect to lux. In the next stage, we derive the general formula for the MC, C^/?''^'' = 

^((?^/9(7p9(7p')Im$^^j((7 + (7',a; + iO)/ijj|c^,q,q'^o- Hereafter, we consider ^'■p} {2q + 2q' , lu x) to simphfy the calculations, 
and divide the final result by four. 

In eg. (|l7|) , the terms with vector potential(s) in j^(g + q',0), given by eq.(^, were dropped. Actually, it is easy 
to see that they make ^ij,v{q + q',oj) independent of q or q' , which lead to C°^f^'^'' = 



C. Vertex Corrections 

The definition of the one particle Green function is 

G(fc,6„) = -i^^drdr'e'^"(--^')(r.c.fe(r)ct,(r')), (18) 

where e„ = vnT{2n+ 1) is the Matsubara frequency for fcrmion. By introducing the self-energy E(fc,e„), the Green 
function for the interacting system is expressed as 

G(fc, e„) ^{e^+^i-el- E(fc, e„))-\ (19) 

where /i is the chemical potential. 

As shown in Fig.^ (i) , the full four-point vertex F is given by the following integral equation: 

F(fc_e+; k+e\k'^e'- k'_e'+) ^ F^(fc_e+; k+e\k'_^e'; fcLe'+) 

+r r'(fc-e+; K^WU"; k'Lt"+)G{k'L,e"+)G{kl,e")T{k'Le"+- kle"\k'^e'- k'_e'+), (20) 

k"e" 

where we put k± — k ± q and — e„ + oj^. F^ is the 'irreducible' four point vertex with respect to the particle-hole 
pair. Here and hereafter, we do not write the spin suffix explicitly, and the suffix of the Matsubara frequencies n and 
A are sometimes dropped for simplicity. 
The current vertex function is given by 

A,(fce+; ke) = vl{k) + r(fce+; fce|fcV; fc'e'+)G(fc', e'+)G(fc', e')vl{k'), (21) 

fc'e' 

as is shown in Fig.|l| (ii). In the same way, 

K{k-e- k+e) = vlik) + F(fc_e; k+e\k'+e'; kLe')G{k'_^,e')G{k'_,e')vl{k'). (22) 
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Aa 



k- k' 



FIG. 2. The definition of {rr^(=')}(fc_;/fc+|fc:i:;fc:^|A;';fc'). 

We also introduce the vertex A^^ which is related to the bare vertex v'^p{k) — dape^ as follows: 
Kp{k^-e; k++e) = vl^{k) + T^^^ik^^; k++\k'^+; k'__) ■ G{k'__)G{k'^+) 

k'e' 

X ( vlpik') + A„(fc^_, fc;_)G(A;V_)A/3(A:V_; fc;+) + Ap{k'__; k'_+)Gik'_+)Ac{k'_+; k'_^+) ) 
+T' ^ {TT'^^^}{k^^;k++\k'^;k^k[;^:k'^_)-G{k'_)AJk'_:k'^^^^^ (23) 

k' k" e' e" 

where k±± = k ± q ± q' , fc±zp = k ± q ^ q' and fco± = k ± q' , respectively. (Here frequency variables were dropped 
because it would not cause a confusion in that case.) The eq.(p3|) is expressed diagrammatically in Fig.|| (iii). Here, 
r-^^'^) is the irreducible six-point vertex, and the vertex {Fr^^^^} is defined as 

{TT'^^^{k^^-k++\k'^-k'_\k'^^-k'^_) 

= Tj2i 5k,p/T + T{k^^-k++\P++\P—)G{p—)G{p++) )r^(^)(p__;p++|fcV;fc^|fc^V;^o-), (24) 



which is also depicted in Fig. ||. 

In the case of g = 0, Ward identity assures the following relations p7|jl^]: 



Aa(A:e„;/ce„) = w°(fc) + 9aS(A;,e„), 
Aa/3(fce„; ken) = v%{k) + daf3^{k, e„). 



(25) 
(26) 



Now we analyze 'I>|fj(2q + 2q' ,uj\) in terms of the diagrammatic technique. They are expressed in Figs. ^ and IJ. 
Below, we show analytical expressions for some of them: 



(a) = -T^A^(fc++;fc__)G+(fc__)t;«„^(fc)G(fc++), 

k,t 

(b-1) = -T^ A^(fc++; fc__)G+(fc__)A,(fc_„; fc„„)G(fc__)A„;3(fc__; fc++)G(fc++), 

k,t 

-TY,K^{k++- fc__)G+(fc__)A+(fc__; fc+_)G+(fc+_) 



(d) 



fc.e 



xA^(fc+_; fc+_)G(fc+_)A^(fc+_; fc++)G(fc++), 



(27) 
(28) 

(29) 



where = duap^^^ and we write G^{k) = G{k, e+) and A+(fc_; fc+) = Aa{k^e^; fc+e+), respectively. We also give 
expressions for some terms in Fig. ^, which include four, six and eight point vertices: 

(f-l) = -r-'' J2 T"~^\k^^e+-k++e\k'e'-k'e'+\kl+e",k'L_e") 

kk' k" ,ee' e" 

xG{k++)Af,{k++; k^^)G+{k^-) ■ G+{k')A^{k'; k')G{k') 



xG(fc^_)A„0(r_;fc^+)G(fc^+), 
(g-l) = -r2 ^ F(fc__6+;fc++e|fc^.+e';fc^,_e'+) 

kk' ,ee' 

xG(fc_+)A„(fc_+; fc++)G(fc++)A,,(fc++; fc__)G+(fc__) 
X G+ (fc^,_ ) A, (fc^,_ ; fc^,_ )G(fc[,_ ) AMfc[,_ ; fc^,+ )G(fc^,+ ) , 

(h-1) = -T^ {rF^^3^}(fc+-^+;fc++e|fc'e';fc'6'+|fc^,Ve";fc^6") 

kk' k" .ee' e" 



(30) 



(31) 
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FIG. 3. The expressions for $^^2(2g + 2q' ,u)\), which are composed of A, G and 
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FIG. 4. The expressions for ^^^}[2q + 2q' ,u>\), which contains at least one four (or six, eight) point vertex. 
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1 ^ 

Ime =0 



Im(8+a)) =0 



FIG. 5. The analytic regions 1~3 as a function of a complex variable e (or e'). Here we put lmuj\ > 0. From each region, 
g{e„;Lux), A^(e+; e„) and r(e+; e„|e+, ; e„') are analytically continued to be gi{e;uj), J^''(e+; e) and Tij [e"^ ; e\e' ; e'"*") = 1 ~ 3), 
respectively. 

xG(fe++)A^(fc++;fc__)G+(fc__)A+(fc__;fc+_)G+(fc+_) 

xG+(fc')A.(fc'; k')G{k') ■ G(fc(,'_)AMfc^,'_; fc[,V)G(fco+), (32) 
where {rr^^^^} is defined in eq.(|^. 

III. MAGNETOCONDUCTIVITY: WITHOUT VERTEX CORRECTIONS 

In this section, we derive the MC by neglecting all the current vertex corrections in the first place. We will see 
that this approximation gives the consistent result with the relaxation time approximation (RTA). In later sections, 
we derive the exact formula for the MC of order 7^^, by taking all the current vertex corrections into account. 

A. analytic continuation 

In deriving the conductivity, we have to perform the analytic continuation of <I>)j^; (2g + 2g', uj\). For an instructive 
purpose, we briefly review the conductivity cr^iy without the magnetic field. If we neglect all the vertex corrections 
for the current, the current-current correlation without the magnetic filed, <I'^']](w), is simply given by 

^ -T^G(fc,e+)G(fc,6„) •<(fc)^;°(A), (33) 



where e+ = e„ + ujx. After the analytic continuation with respect to Imujx > 0, we get 



th^5i(A,e;w) + (t^l^^t^^^ j 52(A,e;w) -th— ff3(fc,e;w) 



vlik)v'M, (34) 



where e+ = e + lu, and gi{k,e;uj) = G^{k,e+)G^{k,e), g2{k,e;uj) = G^(k,e+)G^{k,e) and g3{k,e;uj) = 
G'^(fc, e+)G^(fc, e), respectively. 

Here G^(e) = G(e + iO) {G^{e) — G(e — iO)) is the retarded (advanced) function. In eq.(Q), the first, second and 
third terms come from the region 1, 2 and 3 in Fig. |^, respectively. The coefficient 2 in front of the fc-summation in 
eq.(|3^) comes from the spin-degeneracy. Hereafter, we take this factor 2 correctly in all the expressions in this article. 

In a Fermi liquid, G^{k, e) is expressed for e ~ and |fc| ^ fcps as 

G"ik,e)^zk{e + ^i~el+irkr\ (35) 

where Zfe = (1 — (3/9e)ReE(fc, e))~^ |e=o is the renormalization factor at the Fermi energy. We also introduced the 
quantities = + ReE(fc,0), 7fc ~ — ImE^(fc,0) > 0, and their 'renormalized quantities' ej. — Zk{tk — f^) + and 
7fc = Zkjk, respectively. 

At sufficiently low temperatures, 'yl ^ T is satisfied in a Fermi liquid because 7^ oc in 3D and 7^ oc T^logT"^ 
in 2D. In this case. 



92{k,e;0)RiTrzkS{€ + fi-el)/jk, (36) 

Zk 



5i(fc,e;0)« (- "-^^-—S- (37) 
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After the integration with respect to e, eq.(]3q) gives the singular term of order 7" ~ r . Whereas eq.(|37|), which 
have no pole in the region Ime > 0, give no singularities after the e-integration. 

In conclusion, the conductivity fi^^ = lim^^^o Ini$^i,(w)/w conies from the term possessing (72-section in eq.(^) as 
follows: 



(-^) ^fe--<(fcK(fc), (38) 



which is proportional to "f~'^. Note that we get the exact expression if we take account of all the current vertex 
correction completely. It is still of order 7"^ because Jk is of order 7*^. In the next subsection, we perform the analytic 
continuation for the MC in a similar way, and we also see that only the term which comes from the region 2 in Fig. 
^ gives the most singular contribution with respect to 'y~^. 



B. Magnetoconductivity 



In this subsection, we perform the analytic continuation and derive the expression for the MC without vertex 
corrections for the current. In this approximation, (i) all the diagrams in Fig. | are dropped because each of them 
contains more than one four (or six, eight) point vertices, and (ii) the three-point vertices A^j-^,) in Fig |3| are replaced 
with the bare ones, v'^(^^)ik). On the other hand, we notice the Ward identities for Aq, and Aq,^ given by eqs. ( |25| ) and 
(p6|), in order to get the gauge-invariant expression for the MR as shown below. 

As explained in the previous subsection, the most singular terms with respect to 7"^ come from the region 2 in 
Fig. H Thus, we replace G(fc,e+) G^{k,e) and Gik.e) C^ik.e), respectively. In this subsection, we take the 
Qp, g'p/ -derivative of <^^'^\2q + 2q' ^ uo\) in order to derive the MC. For this purpose, we use the relation 



A«(fc_e;fc+e) = V;^(fc,e) + 0{q% 
Kf(fc,e) =«^(fc)+9„I](fc,e + i<5) 



(39) 



(40) 



where Va{k,e) = w°(fc) -I- daReY.{k,e) and ^a{k,e) = — 9QlmE(fc,e -I- iS). They are shown in Fig.|| (i). As a result, 
w-linear term of <I>(^)(2g -I- 2q', w + iO) without the vertex corrections is given by 

lim ^^'^\2q + 2q' ,u) + ]ff)/Lj 

-^{vlc.p{k) + {V^p{k)G^{k+M{k++)+vl{k^-)G^{k^-)V^,{k)) 

+ [(V;«(fco-)G^■(fc+-)«!];3(fc+)+^.°^(fc-)G^(fc_+)F„^(fco+)) 
+Kf(fco-)G^(fc+_)t;°(fc+_)G-^(fc+_)T//(fc+) 

+ {V^{k,^)G^{k+^)Vp\k+)G^{k++)vl{k++) + «^(fc__)G«(fc__)Kf (fco-)G^(fc+-)l^/(fc+)) ■ 
"(a,g)^(/3,g')] }• (41) 



To take the g, q'-derivative of eq. ([4l|) , we employ the diagrammatic technique for the preparation of deriving the 
exact formula in the next section. Then, we use the foUowings: 



G«(fc±, e) = G^(fc, e)±qp- V^{k, e)G''{k, ef + 0{q'), 
V^{k±,e)^f,{k,e)±qp-Vp%{k,e) + 0{q'), 



(42) 
(43) 



where V^(fc,e) = d/3V^{k, e) — daVf-f- {k , e) . There, we took the summations of p and p' implicitly. Eq.(^) is shown 
in Fig.| (ii). 

Figure ^ shows the diagrammatic expressions for the q, g'-derivative of (a) and (c-1) in Fig.^. The q-derivative of 
the Green functions increases the number of G's by one, which increases the singularity of the conductivity by 7"^. 
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(i) 



G(k-q) 



G(k+q) 



Aa 



(ii) G(k±q) 



G(k) 
G(k) I 



a + O(q-) 



= G(k) 



G(k)s> 
G(k)>^t 



P xqp + O(q^) 



vp ± iy„ 



FIG. 6. The expansion of (i) Ac,(fc_e; fc+e), and (ii) G{k±) with respect to q. 



(a) 




FIG. 7. The g, g'-derivative of (a) and (c-1) in Fig.^ without vertex corrections. 'R' ('A') represent retarded (advanced) 
Green functions, 'op' on retarded (advanced) Green functions means Vp — ijp {vp + i'jp). Note that the sign of each diagrams 
comes from d/dq and d/dq' . 
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+ (tt" p), (P**P') with sign change 



FIG. 8. 

The diagrams for the MC of order 7"^ which are composed of only G, v and J. (r-l~4) corresponds to A'^' in eq.(|44|). In the 
same way, (s-1,2) corresponds to A^'^\ (t-l~3) corresponds to A^^\ (t-4,5) corresponds to A^'^\ and (t-6) corresponds to A^^\ 
respectively. (r-l)~(s-2) contain neither 7^ nor 70,^3, but (t-l)~(t-6) contain 7^ and/or 7^^. 



In Fig. ^ the coefficient of a diagram and that of a diagram which differ by a ^ p are opposite in sign. This fact 
ensures the gauge-invariance of the expression. 

To derive the formula for MC, we take the g, g'-derivative of the all other terms eg. (^) . Figure || shows the 
obtained terms for the MC of order 7"^, which are the most singular terms as we will see later. Note that the terms 
(t-l~6) become zero in an isotropic system because they contain the momentum derivative of ImE^ along the Fermi 
surface, [see eq.(^6|) in this paper.] However, they will play important roles in rather anisotropic system like high-Tc 
superconductors. 

In Fig. 1^, we have dropped some diagrams which cancel out exactly with their counterparts. We show the 
example of the cancellation in Fig. ^. Also, the less singular terms with respect to 7^^, some examples are shown 
in Fig. |l^ (i), are discarded. Furthermore, we also dropped the terms proportional to |G(fc, e)!"* TmG^"'(A:, e) or to 
|G(/c, e)|'"-ReG^"+^(A;, e): Some of the them are shown in Fig. [l^ (ii). Because they are odd-functions of (e + /j, — e^), 
their contribution for the MC should be vanishingly small after the e and fc-integrations, except that the density of 
states is extraordinarily asymmetric around the chemical potential. 

In this paper, we assume that the magnetic filed is along the z-axis. We can verify that the obtained expression for 
the MC, which is given in Fig. ^, changes its sign with the replacement a p, or /3 ^ p'. Thus, the gauge- invariant 
result of the MC, C^p'^p' = GjJ^^'^w • eapz^pp'z in eq.®, is assured. 

Below, we give the analytical expression for Figjl. 

E (44) 

1,---,5 

A^^\k,t) = \G{k,e)\^ • 2 [{vly - V:^.^Vyy)Jy + 2{VyV^y - VxVyy)Jy^x + '^{VxVxy - VyVxx)Jy,y 
H~ i^yJi/.-xx '^'^x'^yJu.xy ^~ "^x'-^^^yy^ i 

A^'^\k, e) ^ \G{k, e)pRe{G^(e)} • 4 [{vly - v^^Vyy)J^ - {vp^^^^ - 2v,^VyJ^^^y + vlJ^^yy)] ^^^^^ , 
A(^)(fc,e) = -|G(fc,e)|'*Im{G^(fc,e)} ■ A[2{v^Vyy^r, + VyV^xly - v^v^y^y -VyVxy^x)Ju 

+4(?;xWy7y - vl-ix)Ju,x + '^{VxVy-ix - vl-iy)Ju,y] ^^.^^j , 

A^'^\k,e) ^ -|G(fc,e)pIm{G^'^(fc,e)} ■ 4 [2{vxVyy"fx + VyVxxly - VxVxyly - VyVxyjx)Ju 

+ ivllyy - 2VxVy-fxy + vl-fxx)Ju\ , 
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FIG. 9. Examples of the diagrams which vanish identically. 



(i) 



(ii) V, 




+ c.c 




FIG. 10. (i) The diagrams of order 7 ^, which are negligible compared with the terms 0(7 ^). (ii) The diagrams which axe 
odd with respect to {e + fj, — el): They will be vanishingly small after the integrations with respect to e and k. 
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where Vxy = d^Vy = d^yiel + Rel]fc(0)) and '^xy = d^ly = 5:rylml]fe(-i(5). In eq. @, we replaced in eq.Q) 

with Jf^(u){k,e) = Wp(i/)(fc2e), and J^^ 

.xy = Vi,xy{k, e) . Note that we have checked eq.(p^ by taking 



the g, q'-derivative of eq.(41) directly by using MATHEMATICA. In later sections, we study all the vertex corrections 
which are dropped here, and find that each and in eq.(^) have a vertex correction by 722- (see eq.([70|).) 

In eq. (^) , A^-^'' and A^'^'' are expressed diagrammatically as (r-l)~(r-4) and (s-1,2) in Fig. p|, a nd A^^\ A^'^\ and 
are expressed as (t-l~6). The latter contains the momentum derivative of 7(fc,e). In eq.(U4), the factor 2 due 
to the spin degeneracy is taken into account. 

Finally, we perform the e-integration in eq. (ESI) under the assumption that 7* <C T. In this case, the expression 



for the MR, eq. (44), becomes much simpler as follows: In the case of 7^ ^ T, the following replacements are allowed 
in eq. (Q). 



\G\^Jie{G^}^^Zk5{e + iJ,-el)^, 

-3 



'27r 



|G|4lni{G^} ^ nzk5{e + fi ~ el) — , (45) 



\Gf^nz,S{e + f^-el)^^, 

By using these relations (^5|), we see that (r-1) and (s-1) in Fig. ^ are equal, whereas (r-2) and (s-2) cancel out 
completely in the case of 7^ ^ T. After the long but straightforward calculation, we get the following simple result 
for the MC. 

- ■ 3^ E / * (-1) '^*(^ + " - ■ ^(«. X 4, . V {i-fe X 4) . V (:^ 

4wB Jfs 7k dk\\ \ 7k dk\\ \ 7k 
32 ^-'^ I jc I'^kl f ^ f J^t^ 



^ 4i;s /fs ^'^'^ 7k \dk\\ V7k/j ' ^^^^ 

where we used the relation (e^ x t/fc)-V = v^dy — Vydx — Wk\-^- Here fc|| is the momentum along ey = (e^ x Vk), that 
is, fc|| is along the Fermi surface on the xy-plane. We also used the replacement: 

y^ZkSifi^el) ^ [ jl^dSk in 3D. (47) 

V ^FS \vk\ 

On the other hand, in the two-dimensional case. 



J2zkS{pi-el) ^ / ^dfc|| in 2D. 



(48) 



The result ( ^6[ ) is nothing but the result of the relaxation time approximation given by eq. (^) , if we replace 2Tk ^ 7/j ^ . 
In conclusion, the MC derived in this section, Aa^^, is of order 7^^. 

IV. MAGNETOCONDUCTIVITY: WITH FULL VERTEX CORRECTIONS 

In the previous section, we derived the MC by neglecting the vertex corrections as the first step. The obtained 
expression becomes equal to the result of the relaxation time approximation. However, it will not be applicable for 
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k+q k+q 



(i) d/dcj, 











~ d/dkp 











k' k' 
k-q k+q 



(ii) d/dcj, 



J(2, 



k'-q k'+q 
k-q k+q 



(iii) d/dq, 



k"+q 
k -q 



= To (kik') 



,1(3) 



Fp (kik"ik') 



(iv) d/dcj, 



k+q-q' k+q+q' 

k"+q' 



' k"-q' 



k k 



d/dkp 



k' k' 



FIG. 11. g-derivative of several types of the irreducible vertices contained in "I>JS(2g + 2q';ijJx)- (i) r^(fc+; fc+|fc'; fc'), (ii) 
r\k-;k+\k'+;k'_), (iii) r"-^\k-; k+\k' ■ k'\k'+; k")-d'J,G{k"), and (iv) T^^^' (fc+_; fc++ |fc'; fc'|fco+ ; fco-) -a^'G. Only the g-derivative 
of (i) and (iv) can give the 0(7"'') contributions for the MC. Whereas (ii) and (iii) give at most 0(7"^), which will be explained 
in Appendix B. 



correlated electron systems because the vertex corrections are indispensable to satisfy the conservation laws. In this 
section, we derive the exact expression for the MC in the Fermi liquid with arbitrary interaction, up to the most 
divergent contribution with respect to 7"^. 

At first, we take the Qpg^/ -derivative of <i>^^)(2g + 2q\ uj\) which is shown in Figs. || and ^. Now, we should remind 
that all the diagrams are composed of (a) the bare vertices v'^^{k), etc, (b) the two point Green functions 

G(fc, e), and (c) the irreducible four, six and eight vertices F^, F^*^^^ and F^*^^^. (see Fig.|l|.) To obtain the exact result, 
we have to take the qp, g^, -derivative of the every element of (a)-(c) included in $'^^^(2(7 + 2q',LUx) once, by taking 
care of the momentum suffixes. 

The qp(qp/)-derivative of (b) is already shown in Fig. So we consider on the derivative of (c): Figure |ll| shows 
several types of the irreducible vertices which are included in ^'^^^{2q + 2q',uj\). In the case of (i) and (iv), d/dqp 
can be replaced by d/dkp, while it is impossible for (ii) and (iii). In this section, we consider only the Qp-derivative 
of (i) and (iv) and neglect that of (ii) and (iii), because they turn out to give less singular contributions with respect 
to 7^^: We will give the explanation for this fact in Appendix B. 

From now on, we derive {d^ / dqpdq'^,)^^^'' {2q + 2q' ,uo\) by taking care of the definitions of F and A given in Fig. 
|l. This procedure is not easy because we meet hundreds of diagrams in total. For example, in Fig. |l^, we show 'all 
the terms' come from the gg'-derivative of (a) and (d) in Fig. ||, and (g-l~4), (h-3,4) in Fig. ^. The definition of the 
operation ® and in each diagram is given in Fig. As for the q, g'-derivative of F, only (i) and (iv) in Fig.]!!] are 
taken into account. In Fig. |l^, we write its sign and the factor above each diagram. 

In Fig. 13. (a) does not have the momentum derivative of F because we drop the terms with (ii) and (iii) in FigjTl|. 



On the other hand, (d) contains d'p,Ap, (g-l~4) contains 9pF, and (h-3,4) have the six-point vertex F— , respectively. 
Hereafter, we write da^^vij^) = Av.a{k) and dapAp{k) = K^^apik), respectively. Some of the diagrams in Fig. |l^ are 
expressed as 



1st of (a) = Kik)[G ipG+l • F(fc|fc') • [G ^'p,G+]K^^p{k'), (49) 

3rd of (a) - +TY,k^{k)[dp,G ■ dpG+ + {p ^ p')]A.,a/3(fc), (50) 
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k k k k 

(i) r- : tt| I [ = 

k' k' k' k' 

k k 



k k 



rl|>S = (3„+3:)|n~| 



k' k' 



(ii) r^^ 



= a 



I 



k' k' 



k" + 



k k 



k' k' 

FIG. 12. The definitions of r^-{k\k') and T-{k\k'). According to the Ward identity, they are equal to {da + d'c)r^ and 
{da + d'a)T, respectively. 



5th of (a) - -T^ A4fc){G[VG+](3) + [VG](3)G+}A.,a^(fc), (51) 

fe,e 

1st Of (g-1^4) = J2 d^G+] ■ dpT{k\k') ■ [G gG+]A.,p,(A:'), (52) 

kk' 

2nd of (g-l~4) = J2 Kik)[GdpG+]-mk')-[d'AGG+}] 

kk' k" ^ee' e" 

xr(fc'|fc")-[GgG+]A,,p,(A;"), (53) 
1st of (h-3,4) = ^ A4fc)[G a^G+] • r^^(fc|fc') • [G a;G+]A,,,,(A;'), (54) 

where d'^ = d/dk'^ and 9" = d/dk'^. [A Q^B] = A-daB — B-daA, and [• ■ •](3) means to pick up only the terms which 
is proportional to G^. Moreover, we wrote r(fce+; ke\k'e', A:'e'+) = r(fc|fc'), and T—{k\k') in eq.(54) is defined as 

r^^(fc|fc') = T^Y.^i/T + r • GG+)^k;p) ■ t'HpIp') ■ (i/t + gg+ • r)(,.fc,), (55) 
pp' 

where F^— is introduced as 

T^^k\k') = T ^ r^(3)(^|^/|^//) . Q'^G{k"). (56) 

k"t" 

Both T— and V^— are shown by Feynman diagrams in Fig. ^ (i) and (ii), respectively. 

To derive {d^ / dqpdq'^,)'^'^^^ {2q + 2q', ujx), we also have to take the q, q'-derivative for all other diagrams in Figsj^ and 
^. After all, we get hundreds of complicated diagrams for the MC, which is not a practical formula for applications. 
Fortunately, all the diagrams turn out to be collected into small number of simpler diagrams. From now on, we will 
prove that they are equivalent to Fig. |lj together with Fig. ||if Jf_i(u) is replaced with A^(^). The expressions for the 
diagrams in Fig. fsl are given by 



(u-1) = -2r2 J2 A,{k)[G d^G+] ■ d,T{k\k') ■ [G ^'pG+]K.,Ak'). (57) 

kk' ee' 

(v-l) = -2r2 A^(fc)[G i,G+] • dpT{k\k') ■ [d'^G+ ■ d'p,G - d'p,G+ ■ 9^G]A,(fc'). (58) 



fefc' 



In the same way, we get (u-2)= — (u-l)!^^^, (u-3)= — (u-l)|/3^p', and (u-4)= — (u-2)|^^p'. 
The expressions for (r-1) and (s-1) in Fig. after the replacement Jp(,y) ^fiiv) ^-^e given by 

{v-l) =TY,K^{k)[dp,G ■ dpG+]K,^p{k), (59) 
(s-l) = -r^Ap(fc){ [VG+](3)-G+(G+^G) }K{k). (60) 

ke 
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<s-l> <t-2> <u-2> <u-4> <u-4> <v-2> <r-l> <s-l> 



FIG. 13. 

Some examples of ((?^/9gp9gp/)$'^'(2g + 2q' ,u>\) which come from (a) and (d) in Fig.^, and (g-l~4) and (h-3,4) in Fig.^. We 
note that the diagrams which are less singular than 0(7"''), e.g., the first term of (i) and (ii) in Fig.px|, are dropped in this 
figure. Note that the diagrams with the replacement {ct, p) ^ {l3,p') also exist for (g-l~4) and (h-3,4). Other examples are 
shown in Appendix C. 
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®a| j = a| j ^ I I" 

FIG. 14. The definition of the operation © and 



(u-1) (u-2) (u-3) (u-4) 




(+2) (-2) (-2) (+2) (+2) 

FIG. 15. The diagrams for the MC with the q-derivative of four-point vertices. The exact MC of order 7"^ is given by these 
diagrams together with FigM Here, the coefficient ±2 comes from the exchange between [a, p) (/?, p'). 



17 







ep 

t 4:^ 




ep 










+ 

ep' 









'v,p 



(+1) 



ep 



eap 



/% eP(A 



FIG. 16. 



The digram for (u-4) in Fig.^ is expanded by using a Ward identity, eg. (|6l|) . 





FIG. 17. The diagrammatic expression for (i) A^,c«(fc) = dah.v{k) and (ii) h.^,af3{k) = c'a/3A^(A;), respectively. 



The procedure of proof is that (i) we decompose the diagrams with or daf3 in Figs.^ and |l^ by taking several types 
of Ward identities, and (ii) recognize the one to one correspondence between them and (d^ / dqpdq'p, )^^^^ {2q + 2q' , lux) . 

For example, we can show that (u-2) in Fig. |l^, is decomposed as Fig. ^ by using a kind of Ward identity as 
follows: 

dc,r'{k\k') = T'^k\k') - d'^r'{k\k'), 

daT{k\k') = T^k\k') +tY^ r(fc|fc") • d'^{GG} ■ T{k"\k') - d'j:{k\k'), (61) 

k"t" 

which are shown in Fig.|l2|. To prove these identities, we consider the fc-derivative of a function F{k) = 
^j,, r(fc|A;')G(/c'). Then, {d / dka)F{k) is given by taking the /c-derivative of every G (i) on the 'connected line' 
between two fc-points in F{k) supposing that k runs only on the connected line, and (ii) on all the 'closed loops' after 
shifting all the momenta by k at the same time virtually. After all, we see that (d/dka)F(k) is given by replacing 
each one of G(p)'s in F{k) with {d/dpa)G and taking the summation of them. Thus, eq.(0) is proved identically. 
Note that the contribution from (ii) vanishes after all (l^. In the same way, we can derive the Ward identities for 
daA^(k; k) = K^^a and dap^vik] k) = A,y,Q^, which are expressed as Fig.|l^. K^^ap contains the irreducible eight-point 
vertex F^^"*). 

We note that Aj. ^(fc; k) ^ Ay^x{k; k) whereas Aa/sik; k) = Ajj^(fc; k). As for A^, the diagrams with u° on a closed 
loop cancel out identically because of the Ward identity, Aq = da{e\ + S^). On the other hand, K^, ^ d^{e1 + S^) 
because the external frequency uj is finite in the current study. Because of this fact, the diagrams with on a closed 
loop in Ajy do not vanish. 

By using the identity in Fig. ^ (u-4) in Fig. |l^ (with 9q,F or dapT) is decomposed into 42 diagrams which 
are consist of four, six and eight-point vertices (without SqF and dapT). In the same way, we can decompose all 
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other terms in FigM. At the final stage of the proof, we verify the one to one correspondence between them and 
all the terms of {d^dqpdq'p,)<^^'^\2q + 2q' with signs and factors. Note that each term in Fig.^ is included 
in one of the diagrams (r-l)'^(v-5) in Figs. ^ and 15, as denoted under each term. Especially, all the diagrams for 
(u-1), (v-1), (v-4) and (v-5) are already appeared in Fig.^ In Appendix C, we present the group of the terms in 
{d"^ / dqpdq'p,)<^^'^'> {2q + 2q' ,uj\) which gives (u-2~4) and (v-2,3) precisely. In the same way, one can check that the rest 

of terms in {d"^ / dqpdq'^,)^'^'^'' {2q + 2q' ,ujx) give (r-l~4), (s-1,2) and (t-l-G) precisely. 

In conclusion, we have proved that (d^ / dqpdq'p,)^^^\2q + 2q' ^ux) is exactly given by Figs. ^ and [Tsl. In the next 
section, we perform its analytic continuation with respect to uj\ , and derive the exact expression for the MC of order 



V. ANALYTIC CONTINUATION 

In this section, we perform the analytic continuation of {d^ /dqpdq'p,)^^^'' {2q + 2q',uix) obtained in the previous 

sections, and derive the coefficient of the w-linear term, i.e., Cp^'^f' in eg. (p^) . For this purpose, we have to study 
the analytic continuation of the following two types of functions: 

Ll^icox) ^Tj2M^n;Lox) ■ G(e+)G(e„) ■ A,(e„;a;A), (62) 
L%{lox)=T^ A^(e„;c^A)-G(e+)G(e„)-r(e+ e„|e„,,e+)-G(e+)G(e„0-A.(e„';^A), (63) 

where we dropped the momentum suffix for simplicity. Figs . |^ and |l^ correspond to LJJ^ (lox ) and L^^^ i'-^x), respectively. 
We promise that several symbols for the momentum derivative of Green functions are implicit in eqs.(|6^) and (|6^). 

After the analytic continuation with respect to w>,(A > 0), the and e„' -summations are replaced with the 
integrations, depending on the region I'^B in Fig. |^. As a result, we get {x = a,h) 



(64) 



/ ^ l^^^Kl^'H^;^) + (th^ - th^) K;;j^He;u;) - th^i^-(3) (e; c.) 
X^^«(e;c.) = (65) 
K^^He;c^)^jj:\e;u;)g,{e;u;)-T ^ / ^T",, 6'+)g,(e'; a;) (e'; a;), (66) 

k',j = l,2,3 

where we put i,j = 1,2,3, and 7^ j (e+; e|e'; e'^) is given by the analytic continuation of r(e+; e„|e„'; el^,) from the 
region i for e„ and the region j for e„/, respectively. Its explicit expression is given in cq.(12) of Rcf. 0. Moreover, 
Jp\k,e;uj) is defined as 



J«(fc,6;a;)^<(fc)+ ^ / |^7^,(fce+; fce|fcV; fcV+)g,(fc', e'; a;)<(fc'). 



(67) 



fe' j = l,2,3 



Now we put w = in eq.(|67|) because we study the static conductivity a^u) — > 0). In this case, according to the Ward 
identity [|, 

jW(fc,e;0) = z;^(fc,e)-i7,,(fc,e), (68) 
(fc,e;0) = z;^(fc,e)+i7,,(fc,e). (69) 

(2) 

On the other hand, J^{k, e) = (k, e; 0) can be rewritten as 

Jp{k,e)=vp{k,e) + Y, J ^T,2{ke\k'e')\G{k',e')\\^{k',e'), (70) 

k' 

which is derived by using the foUowing relation satisfied for |e| ^ T: 

^°(^)+ E f^'r,f{ke\k'e')g,{k',e';0)vlik')^v,ik,e), (71) 

k', 3 = 1,3'' 
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where T^^^ is the 'irreducible' vertex with respect to (72-section |^. Here we denote %j{k<i] ke\k'e'; k'e') as %j{ke\k' e') 
for simpUcity. 

Based on the above argument, we derive the general expression for the MC up to the most singular contribution 
with respect to 7~^. At first, we consider the terms in Fig. ^, which correspond to L'^^^{uli) in eq.(Q). By doing the 

same argument in the previous sections, the term if^ii^'' in eq.(|6^) gives the most singular term, 0(7^^). Thus, the 
expression for the MC derived from i4r^^'-^\ zicrj^i,, is given by eq.(^^ where J^(fc,e) is replaced with eq.(|70|). This 
expression is still of order 7"^ because (the momentum derivative of) Ji_i(y) is of order 7*^. 

Next, we study the terms in Fig. which are expressed as in eq.(Q). In this case, all the terms except 

i = j = 2 turn out to vanish identically because eq.(|66|) contains the following momentum derivatives of G implicitly: 

[G« d^G^] = 0, (72) 
(a„G«) • (SpG^) ~{a^p)^0. (73) 

In conclusion, we get the expression for MC which comes from Fig.|l^, Aa^^^^, as follows: 

^ol.=B--^^Y. /^(-|)2{|GnimG|J,(..5,-.,5.) + |G|^<)}^^^^^ 

X E / ^r,,{ke\k'e')2[mi^G\Wl^^ + IGI^M/^^)}^^^^^^^, (74) 



where 



Wjp{k,e) = Vx{k,e)J^^y{k,e) - Vy{k,e)J^,x{k,e), 

W^^Hk, e) = ( -«,(fc, e)7^(fc, e) + Vy{k, e)7,(fc, e) ) J^(fc, e), (75) 



where J^{k,e) is given by eq.(^, and J^.x = d^Jp. and J^j^^^y = d^yJfj.- Moreover, d^iel + S^(fc, e)) = Vx{k,e) - 
ijx{k, e). Note that Jx.y ^ Jy,x- In eq. (y^j , the factor 2 due to the spin degeneracy has been included, vb ~ (27r)'^ in 
the cubic lattice. 

Equation (jTJ) can be rewritten by performing the partial integration with respect to k. After the reformation of 
the equation, we get the following simpler form: 

Aa^^ ^-B'^-^^Ej'i (-1) 2{|GnimGK) -2|GPReG^iyf }^^^^^ 

xE / ^r,,{ke\k'e')2{\G\'\lmG\W('^ + \G\^w!;'^} (76) 

The diagrammatic expression is shown in Fig.^ They are of order 7"^ because the factor 7"^ comes from each 
f;2-lines, and is of order 7. 

In conclusion, the general expression for the MC in the Fermi liquid system, which is exact in order 7"'^, is given 
by 

Aaxx = + Aal^, (77) 

where Aa'^^ is given by eq.(|4^) where Jfj,(„) is given by eq. (|70|) . 

Next, we derived the simplified the expression of eq. (^7\j under the condition 7* <C T. In this case, Aa^^ is given by 
eq.(^) where Jkfj,(u) is replaced with Jp{y){k, e = 0) in cq.(^. As for eq. ([76|) , the following replacement are allowed: 
|G|2|ImG| ^ ^lZk5(e + ^i-el)/2-1l, jG^ReG^ ^ -7rzfe(5(e + /i - efc)/47fc, and |G|4 ^ ^Zk5{e + ^i- el)/2-il. Then, 
we can perform the energy integrations in eq.(|7^) easily. 

After the long but straightforward calculation, we obtain the following much compact expression for the MC for 
7* -CT: 



Aa^, ^^B'-^i ^Ldp{k)D4k), (78) 
4ws Jfs \vk\ 7fc 



d,{k)^{Vx{k,0)^^Vyik,0)^^ ('^^^^^^ 



dky ' 'dkxj \ 7(fc,0) 
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+ (a P), (P p' ) with sign cliange 
FIG. 18. 

The expressions for Aa'l^^^, in eq.([76|). They are of order 7~^. 



= W^^i) {k, Oh-\k, 0) + W^''' (k, 0)7-'(fc, 0), 

D^{k)^—y] f ^T2^{kQ\k'e')\Gk'{e')\'d,{k!)+d,Xk) 
""B ^ J 47ri 



1 



dSk 



VB Jfs \Vk' 



^T22{k0\k'e')—d^{k') + d^{k), 



where Jpg dSk represents the two-dimensional integration on the Fermi sm-face. Both Di^i{k) and d^{k) are real because 
is pure imaginary for w = 0. 
On the other hand, the MC within the relaxation time approximation is given by 



RTA 



4wB Jfs Wkllk ^ 



(79) 



dl{k)^[v,{k,Q)-^^-Vy{KQ)-^)j 



7(fc,0) 



(see eq.(^.) By comparing the exact formula ( [ts] ) with eq. (|79|), we see that v^^^^ is replaced with J^^^ due to the 
vertex correction by 722- Moreover, is replaced with due to another vertex 722- Conversely, we get Z\cr^J^ if 
we drop all 722's in the exact expression for the MC, eq. (|78|) . Note that Df^{k) is of order 7"^, as dp(fc) and d^^{k) 
are, because 722 is of order 7^. This fact means that Aaxx given by eq.([78|) is of order 7"'^. 

Before concluding this section, we consider a system with four-fold symmetry around the z-axis (e^ || B). In this 
case, we get the expression 



Aaxx = -B'- — (h ^—d^{k) ■ D^{k)/2, 



4fB Jfs Wkllk 



(80) 



where d±_ = {dx, dy) and D± = {Dx, Dy). Equation ( pOD means that Aaxx is independent of the choice of the x-axis 
on the xy-plane. 



VI. DISCUSSION AND CONCLUSIONS 



A. The Role of the Vertex Corrections 



In this subsection, we analyze the obtained exact formula for MC, eq.([78|), and discuss the role of the vertex 
corrections (VC's) by 722. As shown in eqs. (||), (||) and ([78|), axx, Aaxy and Aaxx have one, two and three 722's, 
respectively. 
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At first, we discuss an isotropic correlated electron system without any Umklapp processes. In this case, Yamada 
and Yosida proved that jj^\k,e;uj) given by eq.(|6^) diverges in the static hmit (w = 0) This fact brings 

the divergence of the static conductivities axx, ^o'xy and Agxx even at finite temperatures, reflecting the momentum 
conservation law. On the other hand, J/j (k, e; u) gives the total current with the 'backflow' in the case of 7^ ^ a; <C T, 
or in the 'zero-sound regime ' 0,0. In this case Jk is finite. According to our formula ([78|), Jk — C-[jv\^, dk = Cijd^, 
and z3k = Cijd^, where Cu is the 'Umklapp coefficient' coming from 0]- Cu is finite for finite lu. This fact 
means that Ap(uj) = up to even in the zero-sound regime in an isotropic system. 

We also discuss the conductivity of an isotropic free electron system with weak disorders, within the self-consistent 
Born approximation. Then, it is known that Cu coming from 722 becomes the 'geometric factor'. In this approxi- 
mation, 7 = Tiimp J2k' l^fe-fe' P^'/Ofe' (0) where Vk is the Fourier transformation of the impurity potential. According to 
the Ward identity, 722 is given by the ladder-type vertex corrections. By solving the Bethe-Salpeter equation, we get 
Cy^ = nimpJ2k' \Vk-k'\'^Trpk'{0){l-cos6k-k')/'2'y- In couclusiou, the conductivities are given by = Cu -(7162/27717), 
Aaxy = Cu • [ne' I ^rr?^'^) and Aoxx — Cy • (77e^/877i'^7^), respectively. As a result, the relations i?H = l/77e and 
Ap = up to S^, which are equal to the results the RTA, hold in the self-consistent Born approximation if we take 
the VC's into account correctly, because of the cancellation of Cu's. 

Finally, we discuss the role of the VC's in anisotropic systems, which will be much important. For example, the 
momentum dependence of Jk in nearly AF Fermi liquids becomes singular singular due to the VC's by 722 |@~@|- This 
fact explains the seemingly non-Fermi liquid behavior of the Hall coefficient in high-Tc cuprates successfully, without 
assuming any non- Fermi liquid ground state [p^ . 

It is natural to expect that the VC's by 722 are also important for the MR. In high-Tc cuprates, the Kohler's rule 
Apj p cx p^-^ is strongly violated in a wide range of temperatures Based on the exact expression derived in this 
paper, we recently studied this long-standing problem by using the spin-fluctuation theory, with including all the VC's 
to keep the conserving laws ||]. We found the approximate relation Apj p oc • p"^ due to the VC's by 722 in the 
presence of AF fluctuations (^af being the AF correlation length). We confirmed that the factor ^^p, which comes 
from the VC's, well explains the violation of the Kohler's rule in high-Tc cuprates. 



B. Concluding Remarks 

In this paper, we derived the general expression for the MC based on the Fermi liquid theory from the Kubo formula. 
We treated the vertex corrections in the exact way by following the Ward identities. The obtained expression is given 
by eq. ([77|) , which is exact up to 0(7^"^). This expression can be simplified as eq. ([zs]) in the case of 7* T . However, 
we should use eq.(|77|) to obtain the reliable result when 7^ <, T is realized like in high-Tc cuprates. In all these 
expressions, we put the factor 2 due to the spin-degeneracy, assuming the paramagnetic state without the magnetic 
field. Fortunately, these obtained expressions do not contain six and eight-point vertices, although the original diagram 
for $(2^ given by Fig. ^ have them. This fact makes the numerical calculation easier. 

It is noteworthy that we can calculate the MC in the conserving way by virtue of the obtained formula (|77| ) or j?^): 
In the purpose, we get at first the proper 'irreducible vertex' through the functional derivative; F-^ = 511/ 5G. Next, we 
perform the analytic continuation and get the irreducible vertex T^o. Finally, we calculate the full vertex correction 
722 through the integral equation like eq. (^) , and input it in eq. (|70|) for J^(fc, e) and in eq.([78|) for £'^(fc). This 
procedure is classified as the 'conserving approximation' developed by Baym and Kadanoff, which is indispensable 
to derive reliable transport phenomena |12| . In future, we can apply the present work to the study of interesting 
magnetotransport phenomena in various kinds of strongly correlated electron systems. 

In the obtained expression, the z-axis is taken to be the direction of the magnetic field B. We obtain the 'transverse' 
MC if we put p = V on the x, y-plane. The transverse MR is given by eq.(|^). In the same way, we get the 'longitudinal' 
MC if we put p = V = z. In this case, the longitudinal MR is given by Ap/ p = ~Aazz/<y1z because of Aazx — 0. 
In both cases, Aa ^ 0{t^) and Ap/ p ^ 0{t^), respectively. In usual three-dimensional systems, transverse MR and 
longitudinal MR are of the same order [|]. In quasi-two dimensional system like high-Tc cuprates, on the other hand, 
the longitudinal MR on the a6-plane in the case of B \\ a, {Apa/ Pa)\\, is much smaller than the transverse one in the 

case of B \\ c, {Apa/ Pa)±- Their ratio will be about 0{t'^/t'^), where tc{ta) is the hopping integral along the c(a)-axis. 

Because we discussed the MC (MR) up to 0{B^), our formula will be valid under the weak magnetic field, 
max/j{ajc/7fc} ^ 1, where Wc = eB/mc is the cyclotron frequency. In general, this condition is satisfied for B <C ITesla 
in good metals at low temperatures. In optimally-doped high-Tc cuprates, however, maxk{^c/7k} ^ 1 will be sat- 
isfied even for B ~ ITesla because maxfc{l/7^} ;< T is expected above Tc. In this article, we did not consider the 
MR caused by the spin-dependence of 7^^^ due to the Zeeman effect, because it is negligible in high-Tc cuprates [ pO| . 
However, this spin-effect may cause the negative MR in the vicinity of the ferromagnetic instability pT[|. 
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Finally, we make short comments on the reliableness of the calculation: We used the MATHEMATICA for the 
check of some parts of the calculation. Moreover, (1) The gauge invariance of the obtained expression for ^Gxx is a 
severe verification of its correctness, because the gauge invariance is violated if we modify only one of the coefficient 
of diagrams for ^a^x- We have to mistake at least four coefficients of diagrams accidentally to recover the gauge 
invariance. (2) The obtained ^Uxx coincides with Acr;^J^ if we drop all the vertex corrections. (3) The exact expression 
for AcTjjy up to O(t^), which was originally obtained in Ref. can be derived much easier by using the technique 
developed in this article. 
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APPENDIX A: THE HAMILTONIAN AND THE CURRENT IN THE MAGNETIC FIELD 

In this appendix, we derive the expression for eqs. (|8|), ( [l^ and (|3|) in §11. Here we study the one-di men sional 
system with only the nearest neighbor hopping for simplicity. It is easy to verify that the obtained results ( A2)-( [A5| ) 
are valid for any dimensional systems with arbitrary range of hoppings. 

In the case of eA <C 1, the Peierls phase factor for ti^j+i, eq.(^, can be expanded as 

exp(±ie(Af°\ + 4°*)/2) = 1 ± + l)e"'' ± '"^^i^"^' + l)^"''' 



p2 

AA'(e'« + l)(e'« +l)e"(9+«) 



(Al) 



By expanding and t*^ in eq. (|6|) in the same way, and by performing the Fourier transformation, we get the following 
result: 

k 

+ e^AaA'pY I - ^ 'Vap,k+(sq+s'q')/2 I ^k+(q+q') /2^k-(q+q') /2- (A2) 
k \ s,s' j 



By expanding u° j,_|_q/2 ^'^^ k±{q±q')l2 ™ eq.(|A^) with respect to qp and q^/, we see that eq.(g) in §11 is exact up 
to 0((j,g') and 0(gg')- 

Next, we consider the current operator in the magnetic field, whose definition is given by eq. (pT|) . We can derive 
that 

j(m) = \\H%,mS^cJ^ 

+ e-^'=(^'"+^'-^)/2c,|„+iC„. - e-'^(^-+^--'/2cLc,„-i} , (A3) 
where the integer m denotes the site index. By performing the Fourier transformation, we get 



d_ 

^ mkk 

_^g-ie(A„+i+A,„)/2g-iA; _ g-ie(A„+yl„_i)/2g-ifc' 



I ^ 4cfe'e'("*=+'''+P)™ |gie(A„.+A„_i)/2gifc _ gie(A„+i+A„)/2gifc' 



'}• (A4) 
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(i) 



d/dK Ip (kik") d/dk; 




(ii) 

m 

(k"lk') 



X" 



d/dki 



(k"lk') 
k' 



(iii) 
^1^ 



(iv) 



A 



p' 



a 

k' 



Tp (kIk') :: 



(kIk') 



r„" (klk"lk'): 



I'p' ^k" 



FIG. 19. The diagrams for the MC of order 7 ^, all of which contain the g-derivative of four-point vertices shown in Fig.pj] 
(ii) or (iii). The terms with the exchanges P ^ p' with negative signs also exist. 



After the straightforward calculation from eq.(A4), we obtain the following expression for j^{p) 



d 



1 "■^'"^ 



d_ 

+ ((a,g) -.(/?, g')) 
d 
dpu 



- 1;° 

P,k+(p+sq)/2 "a0.k-{p+sq)/2 



1 ~ 

k s,s' 



fi,k+{p+sq+s'q')/2 ^a(3,k-(p+sq+s' q') /2 



+ •••, (A5) 

By expanding w° ^±5/2' '"ai3,k±(q±q')/2 ^-nd Vap,k±(p±q±q')/2 '^^ eq.(^) with respect to qp and q'^,, we can check that 
eq.(|l^) in §11 is exact up to 0{q,q') and 0{qq'). 

Finally, we comment that another unquestionable definition for the current operator is j = dHg/dA Based 
on this definition, we can also derive eqs. (|l|) and (|l|) exactly up to 0{q,q') and 0{qq'). 



APPENDIX B: THE MC OF ORDER 7^: THE NEXT SINGULAR CONTRIBUTION 



In this appendix, we study the most divergent terms which come from the q-derivative of the irreducible vertices 
(ii) and (iii) in Fig. |ll|; and Fp^^^ For this purpose, we take the gg'-derivative of ^^'^^{2q + 2q') given in Figs.^ and 

^, and gather all the diagrams which contain F^ or Tp^^\ which have been neglected in the discussion in §IV. This 
procedure produces a huge number of complicated diagrams. Fortunately, we find that they are collected into eight 
diagrams shown in Fig. ^9|, by using the Ward identity. In Fig. |l9|, only the most divergent terms are presented. This 
procedure is similar to that developed in §111. 

For example, the expressions for (i) and (iii) of Fig. are given by, 

(i) =T^Y1 K{k)GG+ ■ d'iTl{k\k") ■ {6k",k'/T + GG+ ■ F) 

kk'k" 

•([G gG+]A^,,,(fc') + {d'p,G ■ d'pG+ - d'^G ■ d'p,G+)Ap{k')), (Bl) 
(iii)=T4 ^ A^{k)GG+ ■T'p<^'Hk\k"W")-d':G-{Sk"M'/T + GG+ -T) 

kk'k"k"' 

i[G d^G+]A^y + {d'p,G ■ d'pG+ ~ d'pG ■ d'p,G+)A^), (B2) 
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+<c.«P> 



(+1) (+1) (+1) (+1) (+1) (-1) (+2) (+1) 




FIG. 20. The part of the diagrams m(d^/dqpdq'p,)$^^^{2q + 2q' ,lo). These diagrams together with the terms with the 
notation '(^-4)' in Figjl^ give (u-4) in Fig.|l5|. 



where and Fp^'^^ are introduced in Fig 



as L^^jj in eq. 



i), the term G+G 



11. After performing the analytic continuation of them in the same way 
qRqA gives the most singular contribution. By using the relation ( p5| ) together 

In conclusion 



we 



with the fact that is of order 7, we can clearly show that eqs.(Bl) and (B2) are of order 7^^ 
can drop Fig. ^safely in usual metals, where 7^ is much smaller than the Fermi energy. 

We can also show that Fig.^ gives the exact 0(7~^)-terms of the MC, and any diagrams which are not included 
in Figs||, |l^ and pj| are at most 0(7"^). We note that the expression in Fig.^ should be gauge-invariant because we 
take all the terms 0(7"^) into account, although it is not written in a gauge-invariant form explicitly. 



APPENDIX C: DERIVATION OF FIGS. | AND ^ 

In order to derive {ff^ / dqpdq'p,)(^^^^ {2q + 2q', u), we have to take the 9, q'-derivative of all the terms Figs. || and |^. 
In this appendix, we present the part of them which give (u-l^A) and (v-1'^5) in Fig. ^ explicitly. Note that (u-1), 
(v-1), (v-4) and (v-5) have already been included in Fig.O. 

At first, we concentrate on (u-4) in Fig.|l^. By using the Ward identity, eq.(|6ll), (u-4) is decomposed into Fig. 
[T6| . There, A^^a and A^^ap are expressed as Fig. |l^, which are also kind of the Ward identities. Thus, (u-4) 
is decomposed into 42 diagrams in total. On the other hand. Fig. ^ show the part of the terms coming from 
[d^ / dqpdq'p,)^^^\2q -f 2q',ui). Then, we can recognize that the diagrams in Fig.^ together with the terms with 
the notation '(u-4)' in Fig. ^ give (u-4) precisely, by checking the one-to-one correspondence between them. In the 
procedure, the symmetry {a,p) ^ {(3,p') should be taken into account carefully. 

In the next stage, we study (u-2), (u-3), (v-2) and (v-3). As for (u-2) and (v-2), they are decomposed into Fig.^by 
using the Ward identity, Eq. (|6l|) . Moreover, A^_c[(;3)(fc) = da{i3)Ap{k) is expressed diagrammatically in Fig.|l^. Thus, 
each (u-2) and (v-2) is decomposed into 9 terms. In the same way, each (u-3) and (v-3) is given by 4 terms. On the 
other hand. Fig. |2^ show the part of the terms coming from {d'^/dqpdq'p,)-^>^'^H2q^2q',uj). Then, we can also see that 
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(-1) (-1) (-2) (-1) (-1) (-1) 




(-1) (-1) (-2) (-2) (-1) (-1) 




FIG. 21. The expanded expressions for (u-2) and (v-2) given in Fig.ha. They are derived by using a Ward identity, eq.(pl| 



ep j 

(h-1,2) :epA 1 (i) : ep| 

d/dkJV 




(-1) 



(-1) 



(h-l,2):epA Ap(i):e|iA ip' 



FIG. 22. The part of the diagrams in {d^/dqpd£i,,)$^^'>{2q + 2q',uj). These dia grams together with the corresponding terms 
in Fig jl| give (u-2), (u-3), (v-2) and (v-3) in Fig.|l5[ 
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(-1) (-1) (-1) (-1) (+2) 




(+2) (+2) (+2) (-4) 



FIG. 23. The rest of the diagrams which come from {d^ /dqpdq'^,)^^^^ (2q + 2q' ,ujx)- After the analytic continuation of these 
terms, we get (r-l~4) (s-1,2) and (t-l~6) in Fig.^. 



Fig.^, together with the diagrams with the notations '(u-2)', '(u-3)', '(v-2}' and '(v-3)' in Fig.|T^, give (u-2), (u-3), 
(v-2) and (v-3) in Figjl^ precisely. Then, we have to take the symmetry {a,p) ^ {P,p') into account carefully. 

In the same way, we can write down all the rest of the terms in (d^ /dqpdq'p,)^^^\2q + 2q',ijj) which are of order 
7~'^. It is rather straightforward to see that they are collected into the 9 terms in Fig.^ precisely. We will not show 
the detail to avoid redundancy. Next, we perform the analytic continuation of Fig. ^ to obtain the expression for the 
MC. By dropping several terms which turn out to be zero, we get the results (r-l'-~^4) (s-1,2) and (t-l~6) in Fig.^ 
which is of order 7~^. 

In conclusion, C^p'^'^p' = {d^ / dqpdq'p,)<S>^^'> {2q + 2q',oj + i0)/8 is exactly given by Figs.| and |l| as for the most 
singular terms with respect to 7^^. Because the MC is given by Aaxx = C^x'^^, the derived expression for the MC 
in the present work, eq.([77|) or eg. ([78|) , is exact up to 0(7"^). 



APPENDIX D: COMMENT ON |COND-MAT/0006028 



In recent works, we find the following relations based on the conserving approximation, Aa^y oc Caf7^^ 01 ^^'^ 
Aa^x oc Caf7^^ @' where ^af is the AF correlation length, and 7 = Iml]fc(— i(5). They are the most divergent terms 
with respect to 7~^. However, in a recent e-preprint [O. Narikiyo, cond-mat / 000602^ , the author strongly claimed 



that the above results are inadequate for explaining the non-Fermi liquid behaviors in high-Tc cuprates. His claims will 
be summarized as follows: (i) The conductivity at finite temperatures given by the conserving approximation remains 
finite even if no Umklapp processes exist, which is a serious unphysical result, (ii) In the conserving approximation, 
the Maki-Thompson process violates the Fermi liquid behavior p (x T^. (iii) The contribution for Aaxy from the 
momentum derivative of 7^ is overlooked in Ref. Q| although it is necessary. 

In this appendix, we show that the statements (i)-(iii) are inappropriate, so the conserving approximation gives reliable 
results for high-Tc cuprates as follows. 

First, we show that the statement (i), which is the main claim of the author, is just false. Actually, we can show 
that the conserving approximation based on the FLEX method produces infinite conductivity in the spherical model 
without Umklapp process: We can prove it straightforwadly by following the Yamada and Yosida's transport theory 
given by Ref. by taking account of both the Maki-Thompson (MT) and the Aslamazov-Larkin (AL) processes. 
On the other hand, only the MT-process becomes crucial if the AF fluctuations are strong like in high-Tc cuprates, as 
shown in Ref. [Q.) In this case, the Bethe-Salpeter eq. (6T4) of Ref. is also satisfied by replacing eq.(6-14) with 

Ao{k, k'- k' + q,k-q)= 7rpfe_,(0)pfe.+,(0)pfc.(0)[(7rT)2 + e^]W,{0), (Dl) 

where Wg(0) is given in eq.(A6) of Ref. jj]. Moreover, ImEfc(— iJ) = i J2k'q ^o(fc, k'; k' + q, k — q) in three dimension. 
As a result, by following the discussion in Ref. ||l^, we see that the FLEX approximation can produce the fact that 
axx = 00 even for T > unless the Umklapp processes exist, as a consequence of the momentum conservation law. 

It is a commonsense that the conserving laws are quite important for the study of transport phenomena; see Ref. 
Ip^ . We found that it is also the case in high-Tc cuprates. 

The statement (ii) is inappropriate: The MT process does not prevent the low temperature Fermi liquid behavior 
p (X because the MT process does not give a singular temperature dependence of Jk at the cold spot. In fact. 



27 



according to the analysis in Ref. Q, Jk Vk/2 at the cold spot due to the MT process is the presence of the strong AF 
fluctuations. Secondly, ImSfc(— i^) oc in 3D is ensured in the FLEX approximation at sufficiently low temperatures 
unless its ground state has a AF order, because the FLEX is classified as the conserving approximation. As a result, 
the MT process in the conserving approximation does not violate the relation p oc in the Fermi liquid regime, (fn 
the numerical results of Ref. |4|, the lowest temperature is too high to see the behavior p oc T^.) 

On the other hand, the Hall coefficient is proportional to the coefficient {dOj /dk\\) at the cold spot. In Refs. 
we find that the MT process produces a seemingly non-Fermi liquid behavior of i?H at higher temperatures where 
^AF is strongly temperature dependent. On the other hand, _Rh should be nearly constant in the Fermi liquid regime 
where ^af becomes nearly constant. 

The statement (iii) is false; actually, Aaxy within the relaxation time approximation is given by 



dk\\ 



FS 



r 9 r 

dkn 



(D2) 



where Ik = Vkj^ ■ If take the vertex corrections into account, Vk is replaced with Jk in eq.(D2). Thus, Aa^y 
seems to contain the momentum derivative of -fk- However, eq.(|D2|) (Ik — Jklk^) is rewritten as 



Agx 



dh 



FS 



dk\\ 



(D3) 



where 0j{k ) = ta.n~^ (Jx / Jy) k = ta.n~^ {lx/ly)k- (see Ref. Thus, the momentum derivative of jk does not apprear 
in eq. ( |D3D . As a result, the calculations of Ru in Refs. based on eq. (D3) are correct. 

In conclusion, our results of the Hall effect and the magnetorestance based on the conserving approximation are 
reliable and consistent with experiments. We find that only the MT-type process plays an important role in the 
transport phenomema in near ly AF Fermi liquid. We note that some vertex corrections are dropped in eq.(3) in 
Ref. [26] of cond-mat / 0006028 , which is for quasiprticle contribution to A(Txx- The correct expression for the MC of 



order 0(7^. ) is given by eq.([77|) in this article, or by eq. (|7^) in the case of 7^ ^ T. 
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